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We present first-principles calculations of the electronic band structure and spin-orbit effects 
in graphene functionalized with methyl molecules in dense and dilute limits. The dense limit is 
represented by a 2 x 2 graphene supercell functionalized with one methyl admolecule. The calculated 
spin-orbit splittings are up to 0.6 meV. The dilute limit is deduced by investigating a large, 7x7, 
supercell with one methyl admolecule. The electronic band structure of this supercell is fitted 
to a symmetry-derived effective Hamiltonian, allowing us to extract specific hopping parameters 
including intrinsic, Rashba, and pseudospin inversion asymmetry spin-orbit terms. These proximity- 
induced spin-orbit parameters have magnitudes of about 1 meV, giant compared to pristine graphene 
whose intrinsic spin-orbit coupling is about 10 (ieV. We find that the origin of this giant local 
enhancement is the sp 3 corrugation and the breaking of local pseudospin inversion symmetry, as 
in the case of hydrogen adatoms. Similarly to hydrogen, also methyl acts as a resonant scatterer, 
with a narrow resonance peak near the charge neutrality point. We also calculate STM-like images 
showing the local charge densities at different energies around methyl on graphene. 

PACS numbers: 72.80.Vp, 71.70.Ej, 73.22.Pr 


I. INTRODUCTION 

Spin-orbit coupling (SOC) effects in graphene func¬ 
tionalized with adatoms and admolecules offer new pos¬ 
sibilities for tailoring and manipulating electron spins, 
potentially leading to new spintronics devices-’’ . It has 
already been demonstrated that a giant enhancement of 
the rather weak intrinsic SOC of the Dirac electrons 
in graphene ! can be achieved by adsorbates, such as 
light “ and heavy “ adatoms. Inducing large SOC in 
graphene is important for studying spin relaxation ’ 
as well as spin transport, in particular the spin-Hall 
effect 19,2 . Here we show that giant SOC can be also 
induced by organic molecules, taking methyl as their rep¬ 
resentative. 

Methyl radical CH3 is the most simple organic molecule. 
Comprising one carbon bound to three hydrogen atoms, 
it forms an important building block for organic com¬ 
pounds. It is a likely contaminant for graphene, espe¬ 
cially in samples prepared by CVD (chemical vapor depo¬ 
sition) during which a H 2 /CH 4 gas mixture is used- . As 
a result, both H and CH 3 impurities could be expected. 
Closer investigations revealed that hydrogen affects the 
thermal stability of CH 3 trapped on graphene !2 ’ 23 , form¬ 
ing clusters at high temperatures. 

There have already been several investigations of methyl 
bonded to graphene, including magnetic and mechan¬ 
ical ' effects. It has been shown by density functional 
calculations, that a large class of organic molecules (in¬ 
cluding methyl) induce a spin-| magnetic moment on 
graphene l , but in general, there is a strong dependence 
of the induced magnetism on the location, distribution 
and coverage of CH 3 on graphene. Another reason to in¬ 
vestigate methyl radicals is their similarity to hydrogen. 
As shown in Ref. (> , an effective Pauling electronegativity 
of 2.28 can be associated to CH 3 , which almost coincides 
with that of hydrogen , 2.20. Thus the bonding be¬ 


havior of these adsorbates should be comparable. The 
important question is, will also the induced spin-orbit 
phenomena be similar? 

In this paper we present first-principles calculations on 
methyl functionalized graphene in two different limits: 
dense and dilute. For the dense limit we present the 
calculated electronic band structure and spin-orbit split¬ 
tings of bands close to the Fermi level. For the dilute 
limit we take a representative 7x7 supercell with a sin¬ 
gle methyl admolecule, where we calculate the electronic 
band structure, fit the bands at the Fermi level to an ef¬ 
fective symmetry-based Hamiltonian and obtain the 
relevant SOC parameters: intrinsic, Rashba, and PIA 
(for pseudospin inversion asymmetry). Further, we in¬ 
vestigate the nature of resonant scattering of a methyl 
group bonded to graphene '. Finally, we provide calcu¬ 
lated scanning tunneling microscopy (STM) images and 
also study the magnetic moment formation by a single 
methyl radical. 

Our main finding is that CH 3 admolecules induce a giant 
local SOC in graphene, by a factor of 100 as compared to 
pristine graphene’s intrinsic SOC of about 10 /ieV ! . In 
the dilute limit, methyl acts as a resonant scatterer, with 
a resonance peak at —8.8 meV below the Dirac point, 
with a full width at half maximum (FWHM) of 4.9 meV. 
Methyl prefers to bind in a configuration in which the 
hydrogen atoms point in the direction of the centers of 
the subjacent graphene honeycombs. What would be the 
energetic cost of a methyl rotation on graphene? We have 
calculated such activation energy and found that a rota¬ 
tion of the methyl group by 60° around its bonding axis 
would require an energy of 0.17 eV. Finally, we find, in 
agreement with previous studies 2 that CH 3 covalently 
bonded to graphene induces a spin-| magnetic moment. 
The paper is organized as follows. In Sec. II we dis¬ 
cuss our calculational methods. Section III presents the 
density functional theory (DFT) results of the electronic 
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properties in the dense limit. Section IV reports on the 
electronic structure and its phenomenological modeling 
for a 7 x 7 supercell with a single methyl admolecule rep¬ 
resenting the dilute limit. Finally, in Sec. V we discuss 
the induced electronic charge density, spin polarization, 
and STM images. 

II. COMPUTATIONAL METHODS 

Our first-principles calculations were carried out 
using the Quantum ESPRESSO suite based on 
density-functional theory ' with plane waves and pseu¬ 
dopotentials : 1 ,,J . We used fully relativistic projec¬ 
tor augmented-wave pseudopotentials with a Perdew- 
Burke-Ernzerhof 34,35 exchange-correlation functional. 
The kinetic energy cutoff for charge density and poten¬ 
tial was 184 Ry, the kinetic energy cutoff for wave func¬ 
tions was 46 Ry, and the convergence threshold for self- 
consistency was 1CV 8 Ry. In general, a fc-point sampling 
of 10 x 10 x 1 was used for self- and non-self-consistent 
calculations, except for the density of states where a 
higher sampling of 20 x 20 x 1 was necessary and for 
the band structure along T-M-K-P we used 80 discrete 
Vpoints. We used a vacuum spacing of 15 A in the z di¬ 
rection to simulate isolated graphene. Spin unpolarized 
ground states were used to study SOC effects. Structural 
relaxations were performed with the Broyden-Fletcher- 
Goldfarb-Shanno quasi-newton algorithm . In Fig. 1(a) 



FIG. 1. (Color online) Structure of methyl functionalized 
graphene, (a) The geometric structure in the vicinity of the 
admolecule, with labels for the methyl group CH 3 , the carbon 
atom that bonds the admolecule CAd, its three nearest, Cnn, 
and six next nearest, Cnnn, neighbors, (b) Unit cell of the 
dense (12.5%) and (c) of the dilute (1%) coverage limit. 

we show the basis of the geometric structure used in our 
calculations, with atomic labels used throughout the pa¬ 
per. The methyl group is labeled as CH 3 , the carbon 
atom that bonds the admolecule as CAd, the nearest car¬ 


bons as Cnn, and the next nearest as Cnnn- In equi¬ 
librium, the hydrogen atoms of CH 3 point towards the 
centers of subjacent graphene hexagons. 


III. DENSE LIMIT 

The dense limit is represented by a 2 x 2 supercell which 
is functionalized with a single methyl group (12.5% cover¬ 
age) [see Fig. 1(b)]. Structural relaxation shows, that for 
the 2x2 supercell the carbon-admolecule CAd - CH3 bond 
length is 1.607 A, the nearest-neighbor CAd - Cnn bond 
length is 1.499 A, and the distance between the next- 
nearest-neighbors Cnnn - Cnnn is 2.481 A [see Fig. 1 (a)]. 
The lattice constant a is 2.479 A, somewhat greater than 
in pristine graphene (2.466 A). Similar to hydrogenated 
graphene, the chemisorption of the methyl group induces 
sp 3 hybridization. The carbon atom CAd, which hosts 
the methyl group, has an out of plane lattice distortion 
A of about 0.355 A. For comparison, in the hydrogenated 
graphene 1, , the lattice distortion is about 0.36 A. 
Bringing into contact an isolated methyl radical and a 
graphene sheet, the methyl radical starts to deform from 
a trigonal planar to a pyramidal configuration. A de¬ 
formation energy of A E = 0.37 eV is needed for such 
structural reconfiguration. In addition, forming the co¬ 
valent bond CAd - CH 3 the graphene carbon atom CAd 
experiences an out of plane distortion A. The bonding 
energy 

E B = -(Ecr+CHa) + ( Eqi + Ecb 3 ), (1) 

is defined as the difference between the ground-state en¬ 
ergy of the methyl functionalized graphene EG r+ cH 3 , an d 
the sum of energies of the deformed (pyramidally restruc¬ 
tured) methyl group Ech 3 and the locally corrugated 
graphene Eqi without the methyl group. Our calculated 
value is Eb = 2.06 eV. Figure 2 shows the first-principles 
computed band structure of the 2x2 supercell. As in 
Refs. 6 ’ 8,28 , a characteristic band appears at the Fermi 
level, which is induced by the methyl impurity. We will 
mainly focus on three bands, which we name conduc¬ 
tion (a), midgap (b), and valence (c) bands, respectively. 
The analysis of ab initio data reveals that states near 
the Fermi level originate mainly from p z orbitals on the 
nearest neighbors Cnn- Bands in the energy windows 
(-20,-15) eV and (—12,-8) eV are mainly from s and 
Px + Py graphene carbon orbitals. They correspond to 
the intact a bands of pristine graphene. In between, from 
— 15 eV to —12 eV, there is a peak in the density of states 
(DOS) which comes from s and p z orbitals on CH 3 and 
CAd, respectively, which is a fingerprint of their covalent 
bonding. The band structure at those energies is less dis¬ 
persive reflecting its molecular states character. The last 
DOS characteristic spans the energy window from —8 eV 
to —3 eV. There contributes mainly CAd carbon with 
its p orbitals which provides an evidence of the internal 
hybridization towards sp 3 structure. Moreover, at those 
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FIG. 2. (Color online) Calculated band structure (left) of 
the methyl functionalized graphene for 2x2 supercell config¬ 
uration, with labels for the conduction (a), midgap (b), and 
valence (c) bands, respectively. The panel at the right shows 
the corresponding density of states for the admolecule and 
relevant carbon atoms. The contributions of different orbitals 
are indicated by the labeled lines. The panel below the band 
structure figure, shows a zoom on the midgap state in the 
energy region from —0.3 eV to 0.3 eV. 



FIG. 3. First-principles computed SOC splittings for 2x2 
methyl functionalized graphene along the T-M-K-F path: con¬ 
duction (a), midgap (b), and valence (c) bands, respectively. 


energies also CH 3 admolecule shows its intrinsic charac¬ 
ter, namely the states participating in bondings among 
the hydrogens’ s and carbon p orbitals. 

Figure 3 shows SOC splittings of the three relevant 
bands, (a)-(c), near the Fermi level. The splitting max¬ 
ima of the valence and conduction bands are nearly equal, 
while the midgap band is split less. However, all three 
are of the order of 0.6 meV. At the time reversal points 
r and M the spin-orbit splittings are zero. 


IV. DILUTE LIMIT 

The dilute limit is represented by a functionalized 7x7 
supercell comprising a single methyl group ( 1 % coverage) 
[see Fig. 1(c)]. Structural relaxation shows that for the 
7x7 supercell the carbon-admolecule CAd - CH 3 bond 
length equals 1.583 A, the nearest-neighbor CAd - Cnn 
bond length equals 1.509 Aand the next-nearest-neighbor 
distance Cnnn - Cnnn equals 2.464 A. The carbon atom 
CAd lias an out of plane lattice distortion A of about 
0.410 A, which is more than in the dense limit. The rea¬ 
son is the more relaxed cr bond network in the dilute limit 
which allows a more ideal sp 3 tetrahedral distortion. In 
both cases, there is a competition between the CAd - CH 3 
and CAd - Cnn bonds in forming the 109.5° tetrahedral 
angle. This competition tends to modify A and also the 
in-plane alignment of the Cnn - Cnnn bonds. In the 
dilute limit the in-plane alignment is geometrically not 
constrained so tightly as in the dense limit case, what 
results in a larger A and shorter Cnnn - Cnnn bond dis¬ 
tance. This behavior is similar, e.g. with hydrogenated 
or fluorinated graphene* ’ . The bonding energy in the 
dilute limit is found to be 2.46 eV, which is larger than 
in the dense case. However, the magnitude is very close 
to the one for hydrogen E B = 2.9 eV (> , supporting the 
fact of a covalent bonding. 

Before we discuss the ab initio results let us introduce 
our minimal tight-binding Hamiltonian. For the descrip¬ 
tion of the orbital part, we employ a nearest neigh¬ 
bor tight-binding Hamiltonian- ’ based on carbon p z 
orbitals since those are mainly contributing to states 
around the Fermi level. Hamiltonian Ttorb consists of 
an on-site energy £ch 3 term for the methyl group, a hy¬ 
bridization T for the hopping between the adsorbate and 
host graphene carbon and the standard nearest neigh¬ 
bor hopping t = 2.6 eV for the remaining carbons in 
the lattice. For simplicity, we model the methyl group 
as a single energy level with one effective p z orbital that 
bonds on top of a carbon atom. The orbital Hamiltonian 
reads as 

^orb £CH 3 

<t a 

-t E + 

-BjSCnn a 

El + c j',( T B i,<?) 

( i,3) 

~^ El E>Lcj, g + C ],cr C t<r); (2) 

(id) 

where (i,j) denotes the summation over the nearest 
neighbors. The operator X\ (X a ) creates (annihilates) 
an electron with spin <7 in the effective p z orbital on 
the methyl group. Similarly, c\ a (c JiCr ) are the creation 
(annihilation) operators for p z orbitals of graphene 
carbon atoms. Specifically, we introduce A\ (A a ) and 
B' icT ( B ia ) as the creation (annihilation) operators for 
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CAd carbon (assuming it is on sublattice A) and its three 
nearest neighbors, Cnn (on sublattice B), respectively. 
The notation and labeling are illustrated in Fig. 4. 


( a ) x O 3 ) £ch 3 



FIG. 4. (Color online) Schematic plot showing notation and 
graphical representation of the minimal tight-binding Hamil¬ 
tonian. (a) Positions and labeling of the relevant atomic sites 
whose p z orbitals enter the model Hamiltonian. Shown are: X 
= CHs, A = CAd, three nearest, Bi, i = 1,2, 3, and six next- 
nearest, Cj, j = 1,..., 6, neighbors, respectively, (b) Sketch of 
the dominant orbital and spin-orbital hoppings near the ad- 
molecule, all carbons in graphene lattice are coupled by the 
nearest-neighbor hopping t = 2.6eV (not shown). 
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FIG. 5. (Color online) Calculated electronic band structure of 
the methyl functionalized graphene in the dilute limit repre¬ 
sented by 7 x 7 supercell. Left panel: First-principles (dotted 
black) band structure along with the tight-binding fit (solid 
blue) for the conduction (a), midgap (b) and valence band (c), 
respectively. Right panel: The corresponding broadened or¬ 
bital resolved density of states for the admolecule and atoms 
in its vicinity, different orbital contributions are indicated by 
the labeled lines. Panel below the band structure figure shows 
a zoom on the midgap band including the tight-binding fit in 
the energy region from —0.02 eV to 0.02 eV. 

Figure 5 shows the DFT calculated spin-unpolarized 
electronic band structure of the fully relaxed 7x7 
supercell (dotted lines) along with the tight-binding 


fits (solid lines). The spectrum in the vicinity of the 
Fermi level shows three characteristic bands, which, in 
analogy with the dense limit, we call conduction (a), 
midgap (b) and, valence (c) bands, respectively. These 
bands can be fitted by two parameters T = 7.6 eV and 
£ch 3 = —0.19 eV that enter the orbital Hamiltonian 
■%orb i n Eq. (2). They were obtained by minimizing the 
least-square differences between the first-principles and 
the tight-binding computed band structures considering 
the three bands around the Fermi level. The shaded 
regions around the K point in Fig. 6 show the k-space 
range employed in the fitting. 

The midgap state becomes more localized, as a conse¬ 
quence of flatter band dispersion when compared to the 
dense limit. This indicates a weaker interaction between 
the supercell periodic images unlike to the case of dense 
functionalization limit. The energy bandwidth over 
which the midgap band extends is here only 10 meV, 
whereas in the dense limit it is 300 meV. The main 
contributions to the three relevant bands (a)-(c) come 
from p z orbitals on CH 3 and the nearest-neighbor carbon 
atoms Cnn! see the orbital resolved density of states 
in the right panel of Fig. 5. This fully acknowledges 
our minimal tight-binding Hamiltonian model which 
implements only p z orbitals. 

To describe SOC effects, we extract the SOC param¬ 
eters from our ab initio data by employing a minimal 
spin-orbit coupling Hamiltonian % so 6 ’ 8 ’ 18 : 


n so = 


iAf 

3^3 

iAf 

3a/3 

2iA K 


} ' (^z)cra c j,C7 + H.C.] 


Cj GCnnn a 


ei Bj aVii (s z ) aC rBj t0 


«*>j» ° 


Y Y X d Aj ) z>aa ,Bj^ +H.C.] 


Bj£ Cnn 

- 2 lA f IA Y Y B iA g X D ij)z,aa'Bp 0 


{(id)) 


iA I 

3^ 


£' E c t 

«*>i» 


a l 'ij(S z ) r j a Cj !a . 


( 3 ) 


Here, symbol s represents the array of Pauli matrices. 
The sign factor v l3 equals —1 (+1) for a (counter-) clock¬ 
wise hopping path connecting next-nearest neighbors. 
Vectors dij and Dij are unit vectors in the xy- plane, 
pointing from site j to i. The last term in Eq. (3) is 
the global intrinsic SOC of graphene with Ai = 12 peV 1 . 
The primed sum therein runs over the sites that are not 
coupled by A(^ nor Af. The fact that the orbital and 
spin-orbital energy scales are different by three orders of 
magnitude, allows us to fit the orbital Hamiltonian % 0 rb, 
Eq. (2), ignoring any SOC contributions. 

The spin-orbit splittings along the high symmetry path 
F-M-K-T within the first Brillouin zone for the three 
bands around the Fermi level are shown in Fig. 6 . The 
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splittings vanish at the time reversal points T and M. The 
maxima of the splittings for the conduction (a), midgap 
(b), and valence band (c) are of the order of 0.1 meV, 
where the largest SOC splitting is experienced by the 
midgap band, which is in contrast to the dense limit 
case. The multiband least-square fits were performed in 
the vicinity of the K point and we extract the following 
SOC-parameters: Af = —0.77 meV, Af = 0.15 meV, 
A r = 1.02 meV and Af IA = —0.69 meV. We observed 
that the main shape of the spin-orbit splitting curves is 
reproduced only by Af IA . Thus, similar to hydrogenated 
graphene' 1 , effects of SOC originate mainly from the 
breaking of local pseudospin inversion symmetry. SOC 
parameters for the methyl functionalized graphene are 
in magnitude comparable with its hydrogenated counter¬ 
part; for comparison, see Table I. 

We stress that only spin-orbit couplings as obtained from 
the fitting of DFT to tight-binding model can be mean¬ 
ingfully compared with the corresponding parameters 
in pristine graphene. The spin-orbit splittings depend 
on the admolecule concentration, and by themselves are 
of little use when compared with the pristine graphene 
or with graphene functionalized by different adsorbates. 
The CH 3 induced SOC parameters, Ai, Ar, Apia, are of 
the order of 1 meV — 100 -times larger than the intrinsic 
SOC parameter Ai ~ 10/xeV characterizing the unper¬ 
turbed graphene. 

Additionally to the 7x7 supercell configuration, we also 


percells are similar, see Table I, confirms that our model 
is robust and reliable for the dilute methyl functional¬ 
ized graphene. Our orbital results are in agreement with 
a 4 x 4 supercell calculation already reported in Ref. . 

Similarity between hydrogenated and methyl function- 
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FIG. 7. (Color online) Left panel: change in DOS A u, Eq. (6), 
for a single impurity limit with parameters T = 7.6 eV and 
£ch 3 = —0.19 eV. The resonance peak appears at E ~ —8.8 
meV with a FWHM ~ 4.9 meV. Right panel: perturbed DOS, 
i/(E) = i/o{E) + t)Av(E), for the admolecule concentration 
g = 0.001% (solid line) and the unperturbed pristine graphene 
DOS (dashed line) near the charge neutrality point. 
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FIG. 6 . (Color online) Calculated spin-orbit splittings along 
the T-M-K-r path for the conduction (a), midgap (b), and 
valence band (c), respectively. First-principles data (dotted) 
are well reproduced by the tight-binding model (solid) with 
Hamiltonian 'Ho rb + 'Hao, Eqs. (2) and (3), using T = 7.6 eV, 
£ch 3 = —0.19 eV, Af = —0.77 meV, Af = 0.15 meV, A r = 
1.02 meV, and Af IA = —0.69 meV. The least square fitting 
was performed in the shaded regions around K. 

calculated a 5 x 5 supercell structure. The first-principles 
data for this case can also be nicely fitted with our tight- 
binding model, however, with slightly modified orbital 
and SOC parameters: T = 7.6 eV, £ch 3 = —0.16 eV, 
Af = —0.39 meV, Af = 0.095 meV, A R = 1.01 meV and 
Af IA = —0.71 meV. The fact that the values for both su- 


alized graphene indicates that the latter should also act 
as a resonant scatterer . To describe the single ad¬ 
molecule limit, we downfold the tight-binding Hamilto¬ 
nian "Herb; Eq. (2), by removing the admolecule p z orbital 
obtaining 8 ’ 18 ’ 38 - 39 


H' loU (E) = Y. a W A ° A ° ( 4 ) 

a 


with 


a(E) = 


T 2 


E - £ch 3 


( 5 ) 


The change in the DOS, Ai /(E), due to a single methyl 
admolecule, is then given by 


A i/(E) = —Im 


a(E) 


d 


t r [1 - a{E)G 0 (E) dE 


G 0 (E) 


( 6 ) 


where Gq(E) is the Green’s function per atom and spin 
for the unperturbed pristine graphene 


G 0 (E) * ^ In 


E 2 


D 2 -E 2 


-iJJ^e(D-\E\) (7) 


with the effective graphene bandwidth D = \J\Z3nt ~ 
6 eV; for details see Refs. v ’ ’ . 

Employing our best-fit orbital tight-binding parameters 
T = 7.6 eV and £ch 3 = —0.19 eV we can investigate 
resonance characteristics of the chemisorbed methyl 


v(E)/[10 atom spin eV] 
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TABLE I. Orbital and spin-orbital tight-binding parameters which fit the band structure for methyl functionalized graphene 
for 5x5 and 7x7 supercells, respectively. Since the different supercell values are comparable, the robustness of the proposed 
tight-binding model is well acknowledged. Comparison of the fitted parameters with hydrogenated and fluorinated graphene 
shows certain similarity between CH 3 and H graphene functionalization. 


X (Adsorbate) 

nxn 

T[eV] 

£x(eV) 

Af(meV) 

Af(meV) 

ApiA(meV) 

A R (meV) 

ch 3 

5x5 

7.6 

-0.16 

-0.39 

0.095 

-0.71 

1.01 


7x7 

7.6 

-0.19 

-0.77 

0.15 

-0.69 

1.02 

H a 

5x5 

7.5 

0.16 

- 0.21 

- 

-0.77 

0.33 

F b 

7x7 

6.1 

-3.3 

_ 

3.2 

7.9 

11.3 


10 x 10 

5.5 

- 2.2 

- 

3.3 

7.3 

11.2 


a taken from Ref. 
b taken from Ref. 


group. Figure 7 shows the change in DOS, A v{E), 
as well as, the resulting perturbed DOS per atom and 
spin, v(E) = iso(E) + r/Au(E), as functions of the Fermi 
energy for the admolecule concentration 77 = 0.001%. 
The quantity i/ 0 (E) = — ^Im{Go(.£?)} = ^ is the DOS, 
per atom and spin, of the unperturbed graphene. We 
clearly see a narrow peak at E ~ — 8.8 meV with a 
FWHM ~ 4.9 meV in Azz(U), indicating that CH 3 acts 
on graphene as a strong resonant scatterer with the 
resonance close to the charge neutrality point. 

Figure 8 shows the dependence of the total energy 



FIG. 8 . (Color online) Dependence of the total energy on the 
rotation angle 1 p of the methyl group. The insets show the 
configuration of the admolecule in the cases ip = 0 ° and 60°. 
Label ”rot” denotes the case where the admolecule was only 
rotated, label ” rot-(-relax” denotes the case of rotation with a 
subsequent relaxation. 

on the angle ip of rotation of the methyl group with 
respect to the reference configuration corresponding to 
the fully relaxed 7x7 supercell configuration used in all 
of our previous SOC calculations. We performed two 
kinds of calculations. First, we rotated the admolecule 


around the CAd - CH 3 bond ( 2 -axis) without structural 
relaxation. Second, we also relaxed the structures at the 
given angle, keeping it fixed to avoid a back rotation 
during the relaxation process. It turns out, that the 
rotation of the methyl group by an angle ip = 60° 
around the 2 -axis requires a maximum energy of about 
0.17 eV; see Fig. 8 . This corresponds to a temperature 
of 2090 K or a frequency of w = 41 THz. Comparing 
the maximum rotational energy with the thermal energy 
at room temperature kpT ss 25 meV, it is unlikely to 
rotate the methyl group at moderate temperatures by 
just thermal excitations. Another possibility would be 
terahertz radiation. However, to couple the terahertz 
radiation to a molecule, a dipole moment oriented in 
a suitable direction needs to be present. As a matter 
of fact, the CH 3 molecule in its pyramidal geometry 
possesses only an effective dipole moment along the 2 
axis. This is due to the different electronegativities 
of hydrogen (2.20) and carbon (2.55), so the rotational 
excitation by terahertz radiation is not likely to be 
observed. 


V. CHARGE DENSITY, SPIN POLARIZATION 
AND STM 

In the left panel of Fig. 9 we show the top view of 
the electronic charge density p(r) = J2 n & \4’n( r )\ 2 that 
is summed over the eigenstates q with energies £„ in 
the energy window e m i n = —0.2 eV and e max = 0.2 eV 
with respect to the Fermi level. The dashed line corre¬ 
sponds to the cross-sectional view displayed at the bot¬ 
tom. One sees a preferential localization of the electronic 
states mainly on the sublattice that is opposite to one 
where the methyl group is chemisorbed. Therefore, the 
midgap band is formed mainly from states of the sublat¬ 
tice that contains Cnn carbon atoms. The charge den¬ 
sity is strongly centered near the impurity, meaning the 
interaction among the different periodic images is neg- 
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FIG. 9. (Color online) Methyl functionalized graphene in 7 x 7 supercell configuration, (a) Top view of the electronic charge 
density. The charge density was obtained by summing the absolute squares of the Kohn-Sham states that lie in the energy 
between —0.2 eV and 0.2 eV with respect to the Fermi level, (b) Top view of the spin-polarization. The spin-polarization was 
obtained by taking the difference between the spin-up and -down densities (see explanation in the text). Dashed lines in the top 
figures show directions for the cross-sectional views displayed at the bottom, the corresponding color values of the isosurfaces, 
in units (A _,i ), are shown beside the cross-sectional views. 


ligible and a 7x7 supercell is sufficient to represent the 
dilute limit. It is worth to mention the pronounced tri¬ 
angular shape of the electronic charge density. Carbon 
atoms that mainly contribute to it are aligned along the 
directions spanned by the hydrogen atoms of the methyl 
group. In general, this triangular anisotropy is universal 
for systems with C^ v symmetry. It is also present in hy¬ 
drogen and fluorine functionalized graphene. 

We also calculated the spin-polarization, A p = p up — 
Pd own 5 for the 7x7 supercell configuration; see right 
panel in Fig. 9. The spin-polarization was obtained by 
taking the difference between spin up, p up , and down, 
Pd own 5 electronic densities. Each particular spin density 
Pup/down (^*) = E n,k l^n.up/downM ? is obtained as a sum 
over the eigenstates with the corresponding spin polar¬ 
ization and energy that is below the Fermi level. The 
dashed line in the top figure corresponds to the cross- 
sectional view displayed at the bottom. Each sublattice 
carries a different spin-polarization and the up contribu¬ 
tions, are stronger than the down ones. This leads to a 
total magnetic moment of this structure of 1 ps , which is 
in agreement with Ref. and also in the line with Lieb’s 
theorem, stating that an imbalance in the sublattice sites, 
leads to a net magnetic moment (see Ref. ). The local 
character of the spin-polarization, preferentially centered 
near the impurity, indicates that this is induced by the 
adsorbate. The contributions from the hydrogen atoms 
of the methyl group are too small to be visible. Apart 
from the charge density, also the spin polarization shares 
a clearly visible triangular shape anisotropy. 


We also performed STM calculations, based on the for¬ 
malism of Tersoff and Hamann 1 , as implemented in 
Quantum-ESPRESSO code. The STM tunneling cur¬ 
rent is expressed as integral of the local density of states 
between the Fermi level Ep and Ep + eU 

I cx ]T |^(r o )| 2 0(E F +eU- E n , k )Q(E n , k - E F ) (8) 

n,k 

with <(>„ being states of the surface in the energy inter¬ 
val Ep and Ep + eU and ro the position of the tip. The 
image provides information about occupied (unoccupied) 
states for negative (positive) biases U that modifies the 
Fermi level. Figure 10 shows calculated STM images for 
different bias voltages U. The fact that we also obtain 
an STM image for zero bias comes from a smearing con¬ 
tribution in energy, which is added to the bias U. There 
are a few other features, which are noticeable. The first 
one is that for small biases (up to ± 300 mV), only states 
of the sublattice opposite to methyl group contribute, as 
we already saw in Fig. 9. Especially, this can be seen by 
looking at the graphene honeycomb, drawn for U = 0 V. 
For higher biases, we see that the other sublattice comes 
into play, but even at ±1 V, some atoms in the vicin¬ 
ity of the admolecule are not as pronounced as others. 
Another feature is that there is almost no difference be¬ 
tween the images (c) and (d), even if the bias is changed 
by 200 mV. This is due to the fact that the DOS in Fig. 5 
has two gaps at energies near the Fermi level and thus 
no additional states are available. Next, we see that the 
images are symmetric with respect to the bias, which is 
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FIG. 10. (Color online) Calculated STM images within Tersoff and Hamann approach . States lying in the energy interval 
between Ef and Ef + eU are taken into account. Subfigures (a)-(i) correspond to biases between —IV and +1 V. The color 
map gives the values of a tunneling current /, as defined in Eq. (8), in arbitrary units. In subfigure (e) one graphene honeycomb 
is indicated as a guide for the eye. 


not surprising, since the DOS shows this behavior, too. 
Finally, what we notice, is again the trigonal anisotropy; 
states sitting in the direction, where the hydrogen atoms 
point out, mainly contribute for small biases. 


VI. SUMMARY 

We have investigated SOC in graphene functionalized 
by the methyl group, a simple admolecule represent¬ 
ing a wide class of organic compounds, analyzing DFT- 
computed electronic band structures in the dense and 
dilute methyl coverage limits. Compared to the pris¬ 
tine graphene, we have found a giant (100 times larger) 
SOC in the methyl functionalized graphene that origi¬ 
nates from a local sp 3 distortion. We have proposed a 


minimal realistic tight-binding model Hamiltonian and 
provided the relevant orbital and spin-orbital parame¬ 
ters that fit ab initio computed band structure in the 
vicinity of the Fermi level. As hydrogen, also the methyl 
group acts near the charge neutrality point of graphene 
as a narrow resonance scatterer. The minimal model 
Hamiltonian including the fitted tight-binding parame¬ 
ters can be used for further investigations of spin relax¬ 
ation and spin transport, including the spin Hall effect 
characteristics that could be measured in graphene func¬ 
tionalized by light organic admolecules. We have also 
analyzed conditions and energy ranges needed for exci¬ 
tations of rotational degrees of freedom of the methyl 
group. Analyzing the calculated local densities of states, 
which simulate STM images, we found that the electronic 
density near the methyl admolecule shows a characteris- 
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tic trigonal anisotropic shape, which could be potentially VII. ACKNOWLEDGMENTS 
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